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Standard puncture initial data have been widely used for numerical binary black hole evolutions 
despite their shortcomings, most notably the inherent lack of gravitational radiation at the initial 
time that is later followed by a burst of spurious radiation. We study the evolution of three alterna- 
tive initial data schemes. Two of the three alternatives are based on post-Newtonian expansions that 
contain realistic gravitational waves. The first scheme is based on a second order post-Newtonian 
expansion in Arnowitt, Deser, and Misner transverse-traceless (ADMTT) gauge that has been re- 
summed to approach standard puncture data at the black holes. The second scheme is based on 
asymptotic matching of the 4-metrics of two tidally perturbed Schwarzschild solutions to a first- 
order post-Newtonian expansion in ADMTT gauge away from the black holes. The final alternative 
is obtained through asymptotic matching of the 4-metrics of two tidally perturbed Schwarzschild 
solutions to a second-order post-Newtonian expansion in harmonic gauge away from the black holes. 
When evolved, the second scheme fails to produce quasicircular orbits (and instead leads to a nearly 
head-on collision). This failure can be traced back to inaccuracies in the extrinsic curvature due 
to low-order matching. More encouraging is that the latter two alternatives lead to quasicircular 
orbits and show gravitational radiation from the onset of the evolution, as well as a reduction of 
spurious radiation. Current deficiencies compared to standard punctures data include more eccen- 
tric trajectories during the inspiral and larger constraint violations, since the alternative data sets 
are only approximate solutions of Einstein's equations. The eccentricity problem can be ameliorated 
by adjusting the initial momentum parameters. 
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PACS numbers: 04.25.dg, 04.25.Nx, 04.30.Db, 04.30.Tv, 



I. INTRODUCTION 



One of the best candidates for observation of gravita- 
tional radiation is the coalescence of two black holes. The 
successful detection of these gravitational waves by any 
ground or space-based interferometer (e.g. LIGO [l|, 
VIRGO HQ, GEO600 @, eLISA/NGO ||) will be con- 
firmed by matched filtering of the observed signal against 
an extensive compilation of waveforms produced numeri- 
cally from the binary black hole parameter space. There- 
fore, it is important that this collection of constructed 
waveforms represent the most physically viable scenarios 
of binary black hole evolution. 

Full numerical evolution of the Einstein equations 
made possible since 2005 0-Q is the method of choice 
for producing waveforms than span part of the inspiral, 
and the subsequent merger and ringdown phases. The 
waveform templates are dependent on the physical accu- 
racy of the data being used. Currently, the most widely 
used approach for initial data construction is known as 
the "puncture" method. Developed by [Toj], this method 
is able to set up data for binary black holes with generic 
momenta and spins. These initial data are used almost 
exclusively by all groups that evolve black holes using 
the BSSNOK formulation [ll[ together with the moving 
punctures approach Standard puncture initial data 

have the disadvantage that they contain no realistic grav- 
itational waves, because a conformally flat metric is used 
in this method. When evolved, these data lead to a grav- 
itational wave signal that is zero for some time followed 
by a spurious burst. Only after this burst do we get an 
astrophysically realistic chirp signal, thereby making the 



prescription inherently unphysical. Other methods have 
begun to emerge whereby post-Newtonian (PN) approx- 
imations are used, either as the sole contributor [Pa - flij 
or in conjunction with analytical solutions near the black 
hole [15T - ll7l ]. to develop initial data. These methods al- 
low for increased astrophysical accuracy by including, 
to leading PN order, gravitational waves on the space- 
time from the onset of the simulation. In this paper, 
we investigate three methods of initial data construction 
and use them for full numerical evolution of binary black 
holes. These methods can be considered alternatives to 
the standard puncture technique for creating templates. 
We find that evolutions of two of these methods lead 
to improved gravitational waves. There is a chirp signal 
right from the start, and the (still present) burst of spuri- 
ous radiation has a smaller amplitude than for evolutions 
starting from standard puncture data. The burst occurs 
when the built-in post-Newtonian waveform transitions 
to a fully numerically generated signal. 

Throughout this paper, we will use units where G = 
c = 1. Latin indices such as i run from 1 to 3 and de- 
note spatial indices, while Greek indices such as jj, run 
from to 3 and denote spacetime indices. The paper is 
organized as follows: In Section UH we briefly recall the 
3+1 formulation for numerical relativity. In Section IlIIl 
we summarize the standard puncture technique for creat- 
ing initial data and discuss weaknesses in this approach. 
In Section IIV1 we describe our first alternative scheme; 
a post-Newtonian expansion of the 3-metric and extrin- 
sic curvature in Arnowitt, Deser, and Misner transverse- 
traceless (ADMTT) gauge 
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resummed to approach 
standard puncture data at the black holes. In Section fVl 
we describe the technique of asymptotic matching the 
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4-metrics of tidally perturbed Schwarzschild black hole 
solutions to post-Newtonian expansions away from the 
black holes. Also within this section, we describe two 
more initial data schemes: asymptotic matching with 
a low-order post-Newtonian expansion of the 4-metric 
in ADMTT gauge away from the holes, and asymp- 
totic matching of the 4-metrics with a second-order post- 
Newtonian expansion in harmonic gauge away from the 
black holes. In Section IVI1 we describe additional ad- 
justments required for numerical evolution. In particu- 
lar, we describe an algorithm we use to fill black hole 
interiors for initial data that contain singularities. We 
then show comparisons of evolutions starting from dif- 
ferent initial data schemes. Section IVIII summarizes our 
results and includes discussion of future improvements to 
certain schemes. 



II. 3+1 FORMULATION 



Brandt & Brugmann [T(| using Bowen-York extrinsic 
curvature. This method uses a conformal transverse- 
traceless decomposition of the extrinsic curvature [13, HH 
to simplify the constraint equations. A conformally flat 
three-metric, an, is defined as 



(7) 



where ip is the conformal factor. The extrinsic curvature 
is given by 



K ij = ^~ W A% 
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where 



A by = VTiPA n A + V\n\ - {gV - n l A n J A )p k A n A k] 
A=i Zr A 

(t lkl SAkn A in A + e jkl s Ak n M n A ^j , (9) 



For numerical evolutions of general relativity, it is use- 
ful to split the four-dimensional spacctime described in 
terms of the metric g^ v into three-dimensional spatial 
hypersurfaces and time. The Arnowitt-Deser-Misner [l9[ 
decomposition of the Einstein equations defines the line 
element as 

ds 2 = -a 2 dt 2 + gij (dx l + /3 l )(dx 3 + /3 j ), (I) 

where a is the lapse function, f3 l is the shift vector, and 
gij is the 3-metric. The extrinsic curvature is defined as 



K ij = ~2^( dt 9v ~ C P9i 3 )- 



(2) 



Under the Arnowitt-Deser-Misner decomposition, the 
Einstein equations yield evolution equations of the form 

d t gij = -2aK l:j + Cpgij (3) 
d t Kij = (.(/,',, 2/x,,A"i l),D,n • C.l\,, 

SnSij + 4ir gij (S - p) (4) 

and the Hamiltonian and momentum constraint equa- 
tions 



R - K l3 K 13 



+ K 2 = 167T/9 

g»K) = >-f. 



(5) 
(G) 



Here, Rij and R are the Ricci tensor and scalar of the 3- 
metric g^, and Di is the derivative operator compatible 
with g^. In binary black hole simulations, these equa- 
tions simplify due to the vanishing of the source terms p, 
j\ Sij, and S = g l] Sij. 



III. STANDARD PUNCTURE INITIAL DATA 

The most common approach to producing initial data 
used currently is the standard puncture prescription of 



is of Bowen-York form and where n a A is the radial normal 
vector, p A is the linear momentum, and s A is the spin 
of black hole A. This Kij is chosen so that it already 
satisfies the momentum constraint for arbitrary spin and 
momentum. 

The exact solution for a black hole pair with zero mo- 
menta and spins is the Brill-Lindquist solution [22[ 
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where wia is the "bare" puncture mass of black hole A 
and r A = \x — x A \ is the distance from the black hole 
at position x A . As one can see, this conformal factor 
diverges at the puncture location x A . Such divergences 
also occur for non-zero masses and spins. To avoid nu- 
merical problems with these divergences when solving the 
Hamiltonian constraint equation, one usually splits the 
conformal factor [l(| 



t/> = IpBL + U 
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into a regular piece u and the divergent piece ipBL- By 
doing this, the Hamiltonian constraint is reduced to a 
single elliptic equation 
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(12) 

without any divergences, that can be easily evaluated by 
elliptic solvers such as in [23J. Notice that Ag Y = 
for zero momentum and spin, so that u — 0. Further- 
more, if one of the masses (e.g. 777.2 ) is zero we obtain the 
Schwarzschild solution in isotropic coordinates. 

The fact that data are built using the assumption of 
conformal flatness is an inherent problem. It is unrealis- 
tic to demand conformal flatness in the context of astro- 
physical constructions, because it implies that there are 
no gravitational waves on the initial time slice. 
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Often, gravitational waveforms are described in terms 
of the Newman- Penrose scalar, ^4, extracted at some dis- 
tance from the sources. Usually this ^4 is decomposed 
into individual £ and m modes by projecting onto spheri- 
cal harmonics Y^ of spin weight -2 [24]]. Figure [T] shows 
the real part of the dominant t — 2, m = 2 waveform 
mode of ip4 over time extracted at a distance of 90 M. 
We observe a lack of any gravitational radiation until 
~ 80M; this delay in the waveform signal at arbitrary 
extraction radii is common of standard punctures. Also, 
we see the beginning of an outward traveling wave profile 
only after about ~ 15071/ and differing significantly from 
the initial burst; the waves up until this outward pro- 
file region have been commonly referred to as spurious 
radiation. 

If such a simulation is used to make gravitational wave 
templates, the early part of the wave signal has to be 
discarded up to the time when the spurious radiation 
has become sufficiently small. How much exactly one 
has to discard depends on the desired accuracy, but in 
our experience it typically amounts to between 10% and 
20% of the total simulation time. This kind of waste will 
only increase with future code upgrades that will allow 
us to extract the waves at larger distances, since then 
we have to wait longer for the spurious burst to even 
reach the extraction distance. These problems are the 
reason why we seek alternative constructions of binary 
black hole initial data. In [25| it was shown that spurious 
radiation can be reduced by using a 3-metric that is not 
conformally flat. 
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FIG. 1: The real part of the dominant £ — 2, m = 2 waveform 
mode of ip4 over time extracted at a distance of r — 90M. The 
lack of a gravitational wave signature for the initial phase of 
the evolution is a direct consequence of prescribing conformal 
flatness. The burst behavior (visible as the waveform tran- 
sitions from zero to a regular chirp signal) is indicative of a 
so-called spurious radiation signal. 



IV. INITIAL DATA BASED ON 
POST-NEWTONIAN EXPANSION IN ADMTT 
GAUGE 

One alternative to standard puncture data has been 
introduced by Tichy et al. in [12} with further devel- 
opments in fl3l . [Til ] . This method uses post-Newtonian 
expansions in ADMTT gauge |l8[ for the 3-metric and 
extrinsic curvature on the initial time slice. While PN 
theory is strictly valid only away from regions of strong 
gravity such as black holes, it is possible to rewrite the 3- 
metric and extrinsic curvature in ADMTT gauge in such 
a way that they approach the puncture form as in Eq. (|7| 
and Eq. (|9|). In this way one can obtain globally valid 
data provided the initial black hole separation, ri2, is 
within the PN regime. This is seen as a first step away 
from the standard puncture approach and a step towards 
constructing astrophysically realistic initial data. Using 
ADMTT gauge, the expressions for the 3-metric and the 
extrinsic curvature are taken up to 0(v/c) 5 and 0(i?/c) 4 
respectively, where v ~ \J M / r±2 with M = mi + mi and 
c is the speed of light. A formal expansion parameter, 
e ~ (v/c), is used below to distinguish terms of different 
PN order. The 3-metric in ADMTT gauge is given as 
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0(e 5 ). 



(13) 



The full expression for hjj has been computed by Kelly 
et al. [H[ . It denotes the gravitational wave contribution 
to the 3-metric. The conformal factor 
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0(e 6 ), (14) 



is very similar to that of standard punctures in Eqs. (|10[) 
and (fTTj) . The only difference is that it contains an addi- 
tional term at order 0(v / 'c) 4 and no u piece as in Eq. (|11[) 
since the constraints are satisfied only approximately. It 
must be noted that the PN expressions in Eqs. (fTB"]) and 
(|15p are not pure Taylor expansions in (v/c) since terms 
such as ip 4 hi Eq. (|13p are not expanded in powers of e. 
This amounts to adding specific higher-order terms for 
the purpose of creating a metric that is similar to the 
standard puncture expression. This selective inclusion of 
terms does not improve the accuracy of the PN expan- 
sion, however it is performed to ensure the presence of 
black hole apparent horizons in the initial data [l2[ • The 
extrinsic curvature is 
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, TT 
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(15) 

where the leading order term contains the Bowen-York 
extrinsic curvature as in the case of standard punctures. 
The additional terms at order 0(v/c) 5 can be found 
in [l2| • The trace of the extrinsic curvature can be shown 
to vanish up to 0(v/c) . 
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V. ASYMPTOTIC MATCHING OF 4-METRICS: 
AN OVERVIEW 

Another alternative to the standard puncture method 
involves a post-Newtonian expansion of the 4-metric 
away from the black holes. This expansion is then asymp- 
totically matched to tidally perturbed black hole metrics 
used for, and in close proximity to, the black holes. This 
approach is believed to produce more astrophysically re- 
alistic initial data based on several reasons; the method 
is built on physical approximations and not on the poten- 
tially unphysical ansatz of assuming the punctures form 
close to each black hole, and the level of analytic control 
over the physical approximations allows for the method, 
in principle, to accommodate ever higher-order expan- 
sions if desired. 

Matching requires dividing the spacetime into 4 zones; 
2 inner zones, a near zone, and a far zone. The in- 
ner zones are considered as regions where perturbative 
Schwarzschild solutions are valid. The near zone is the 
region where PN theory will hold (provided (v/c) -C 1). 
The far zone is the region of spacetime where retarda- 
tion effects will influence the system. Within each zone, 
4-metric approximations are constructed in some coordi- 
nate system with corresponding parameters (mass, mo- 
mentum, etc.). Matching these different approximations 
becomes possible because their regions of validity over- 
lap in certain regions called buffer zones. Figure [2] shows 
a schematic sketch of the different regions. Asymptotic 
matching [26 29] involves comparing two asymptotic so- 
lutions inside of a 4- volume. This returns a map between 
the local coordinates and parameters of each of the differ- 
ent regions, and forces both solutions to be asymptotic to 
each other within the buffer zone. It is inside the buffer 
zone where these solutions can be merged into a smooth 
global metric by way of a transition function. We refer 
the reader to |15l - fl7| for our choice of transition func- 
tions. 

For a coordinate separation of r\i the gravitational 
wavelength is given by 

A G w = W— • ( 16 ) 
V m 

A diagram of the regions is shown in Figure [5J and a 
summary of the regions' influences is seen in Table 1. 
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TABLE I: Location of inner and outer boundaries of each 
spacetime zone. 




FIG. 2: This figure illustrates the separation of the space- 
time into specific zones: the inner zones, near zone and far 
zone. In each zone one can use a specific approximation. The 
corresponding buffer zones are the regions where two zones 
overlap. 



A. Matching an ADMTT-Gauge Post-Newtonian 
Metric to a Black Hole Perturbative Metric 

One of the first implementations that has used the 
above asymptotic matching approach is described in fioj j . 
In this initial data scheme, a PN 4-metric in ADMTT co- 
ordinates [3, HO, Hl| is used in the near zone. Since 
in this approach the Post-Newtonian metric is accu- 
rate only up to terms of order 0(v/c) 2 , there are no 
retardation effects and the near zone results are also 
valid in the far zone. In the inner zones, a perturbed 
Schwarzschild metric in isotropic coordinates is used. A 
globally valid 4-metric is then obtained by asymptotic 
matching. The matching procedure is simplified by the 
fact that ADMTT and isotropic coordinates are similar. 
The matched result has errors of 0(v/c) 4 away from the 
black holes and errors of order 0(v/c) 2 0(rA/ri2) 3 near 
each black hole. The PN part of the matched results that 
is valid away from the black holes has the same 3-metric 
and extrinsic curvature as in Eqs. (fT3| and (fT5|). but 
with all terms of order e 4 and higher dropped. The per- 
turbed black hole solutions that are valid near each hole 
can be found in [1 61 ] - Both metrics are matched in the 
buffer zones and smoothed through the transition func- 
tion mentioned earlier. 



B. Matching a Harmonic- Gauge Post-Newtonian 
Metric to a Black Hole Perturbative Metric at 
Higher Order 

The second matching approach we consider here is 
from work by Johnson-McDaniel et al. [l7| ■ It applies 
higher-order asymptotic matching of the 4-metric solu- 
tions in the buffer zones. For this approach the inner 
zones approximate the 4-metric in a quasi-Cartesian form 
of Cook-Scheel harmonic coordinates [32j |. the near zone 
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now approximates the 4-metric in harmonic coordinates, 
and the far zone region is now considered through ap- 
proximating the 4-metric in harmonic coordinates. The 
reader is referred to [17| for the complete construction 
of the 4-metric within the near zone regime. The in- 
ner zone solution is expressed to higher-order than what 
was used in the previous matching approach in ADMTT 
gauge and is found in its full form in [17| . The inclusion 
of a far zone contribution introduces an additional buffer 
zone, wherein another transition function will smoothly 
join the near and far zone approximations in the region 
where the errors from both approximations are compa- 
rable. 

The main difference between this matching approach 
and the one prior is that it is now carried out to higher or- 
der; in both the PN expansion of the 4-metric away from 
the black holes and also the perturbation expansion of 
the inner zone black hole solutions. Also, the PN expan- 
sions are now in harmonic gauge. The matched result has 
errors of 0(v/c) 5 away from the black holes and errors of 
order 0(v / c) 5 0(r / V12) 3 near each black hole. This result 
is then further augmented by additional (but formally 
unmatched) PN terms away from the black holes which 
include gravitational radiation, so that away from the 
black holes one has only an error of 0(v/c) e . The choice 
of harmonic coordinates is one of convenience; matching 
has already been established between the near and far 
zones through [33j which simplifies the inclusion of the 
far zone region. This high-order matching yields more 
accurate results than the previous lower-order matching 
in ADMTT gauge. With the inclusion of initial gravita- 
tional radiation, we can expect improvements in the data 
set over the standard punctures approach. 

In [34l| the tidal deformations present in the inner zones 
of [ill (but not the near and far zones) are added to a 
superposition of two spinning black holes in Kerr-Schild 
coordinates. Evolutions display a reduction in the (2,0) 
mode of the outgoing gravitational waves. The expec- 
tation is that a full inclusion of outer zone phenomena, 
as in our work, will also reduce the (2,2) mode spurious 
signal. 

VI. NUMERICAL EVOLUTION RESULTS 

Evolutions from all initial data were done using the 
BAM code [24l l35l - [37l| . The gravitational fields are evolved 
using the BSSNOK formalism [n], |38|, |39[ in the variation 
known as the "moving punctures" method [8|, [9j . Thus 
the 3-metric gij is written as 

9ij = X^% 3 ( 17 ) 

where the conformal metric 7^ has unit determinant. In 
addition, the extra variable 

fi = f^ki-^ (lg) 

is introduced where 7 y is the inverse of the conformal 
metric. Furthermore, the extrinsic curvature is split into 



its trace-free part, Ay, and its trace, K, which is given 
by 

Kii = X _1 Ua + y 7 yj • (19) 

The particulars of our BSSNOK implementation can be 
found in Q. 

All our evolutions use equal-mass non-spinning black 
hole binaries at initial separation of 10M. Also, all 
sets are tested through evolutions using a resolution of 
~ 3M/224 near the black holes. Furthermore, all sets 
use the same expressions for a and j3 l at the initial time; 

a = tPbl' 2 , (20) 
& = 0. (21) 

To evolve the lapse and shift we use the "1+log" slicing 
condition 40] modified with the addition of an advec- 
tion term as used in 0, |4l|, and the "gamma freezing 
condition" [42| modified and used in [2J, |4l|, where it 
is labeled as the "shifting-shift case" in the former and 
the "000" shift choice in the latter. For the initial data 
obtained by matching or from PN approaches, the con- 
straints are only approximately satisfied. In the case of 
standard punctures the constraints are solved numeri- 
cally using the approach in [23]. 



A. Filling Black Hole Interiors 

In the case of initial data obtained by asymptotic 
matching, each black hole is described by a perturbed 
black hole solution. This black hole solution can con- 
tain physical singularities inside the event horizons that 
have to be dealt with before a numerical evolution can 
be attempted. We have developed an algorithm (called 
BHf iller in our code) that fills a specified region inside 
the black hole apparent horizon with smooth data. The 
justification for this procedure is that the physics outside 
a black hole remains unaffected if we only change the 
inside of the horizon. Our algorithm modifies the BSS- 
NOK variables at the initial time as follows: For each 
black hole (centered at the point xa) we pick a sphere of 
radius rfm contained inside the black hole horizon. The 
radius rfm is chosen such that this sphere contains all 
singular points. To set valid data at each point x inside 
this sphere we use a weighted average of linear extrapo- 
lation and a value that corresponds to standard puncture 
data. Let us define r — \x — x*a\, n — {x — XA)/r and de- 
note a particular BSSNOK variable component at point 
x by u(x). If r < rfm we set 

u{x) = u r T + u (l - T), (22) 

where 

u r = u(rfuih) + [d r u(rfiun)](r - r f m), (23) 
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and the weight is given by 



r = J<, 



tanh 
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We apply this method of filling the region inside radius 
Tfui to all BSSNOK variables except T z , which we simply 
recompute from the filled jij using Eq. (|T8| . We only 
need to use this algorithm for the case of asymptotically 
matched data sets, as all other data sets are free of any 
physical singularities. 

Note that the constraints are generically not satisfied 
in the filled region. We have tested our filling algorithm 
by evolving standard puncture data with filling, applied 
at the initial time or some predetermined later time, and 
without filling and compared various quantities. Gauge 
invariant quantities, such as the gravitational wave am- 
plitude as a function of gravitational wave frequency, are 
unchanged when we compare with and without black hole 
filling. We have also checked that in our evolutions no vis- 
ible constraint violations are emitted by the black holes 
filled with this algorithm. The latter is expected since 
the BSSNOK system together with the gauge conditions 
used here has been shown [43[ to lead to causal constraint 
propagation, so that any constraint violations introduced 
inside the black holes cannot affect the exterior space- 
time. Notice, however, that the BSSNOK system has 
superluminal gauge modes (43[. Thus gauge dependent 
quantities such as the lapse are somewhat different, even 
outside the black holes. Hence the black hole trajecto- 
ries (while qualitatively the same) are slightly different 
as well. 



B. Shortcomings of Low-Order Asymptotic 
Matching 

The real part of the dominant I = 2, m = 2 wave- 
form mode of '04 illustrates whether a gravitational wave 
signal has been implemented starting from the t = 
slice. Figure [3] shows the waveforms, extracted at ra- 
dius r — 90M, for a simulation starting from standard 
puncture data and that produced from the first data 
set; asymptotic matching of a low-order PN expansion 
in ADMTT-Gauge to a perturbative black hole solution 
(hereafter called ADMTT- match). As with standard 
punctures, ADMTT-match fails to incorporate gravita- 
tional radiation at t = 0. This is not surprising as the 
3-metric is still conformally flat. 

Another undesired result is that the binary merges 
within a time of 250M, which is much too fast. This be- 
havior indicates an accelerated evolution timescale that is 
not expected for a quasicircular inspiral. Figure 2] shows 
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FIG. 3: The (2,2) mode of 04 extracted at a radius of r = 
90M for standard punctures and ADMTT-match initial data 
sets. Both waveforms show no initial wave signature at t = 0, 
due to their reliance on conformal flatness for the 3-metric 
expressions. It can also be seen that ADMTT-match merges 
much faster (within 250M). This is evidence that the orbits 
were not quasicircular at all. 



the orbit of the binary system in the ADMTT-match 
approach. We see that ADMTT-match leads to nearly 
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FIG. 4: Trajectories of both black holes in the ADMTT- 
match approach. A nearly head-on collision between the two 
holes indicates there is insufficient tangential momenta to sat- 
isfy a PN quasicircular orbit starting from t = 0. 

head-on trajectories for the black holes. The individual 
black holes do not have enough momentum to complete 
even one orbit. One can track this problem back to the 
order of asymptotic matching used in the construction 
of these initial data. The terms in the extrinsic curva- 
ture that are associated with momentum are computed 
from the 4-metric perturbation expansions used near the 
black holes. However, as discussed in the Appendix, the 
matching procedure used in [l6| is of an order that is 
too low to correctly obtain the the 0(v/c) 3 piece of the 
extrinsic curvature that would lead to an orbiting black 
hole. Thus we see that this data set is merely a proof 
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of concept. To obtain useful initial data, we need to use 
higher order matching. 



C. Analysis of Higher-Order Methods 

In this subsection, we compare evolutions starting from 
data obtained by the higher order matching approach dis- 
cussed in Sec. IV Bl and data from the pure PN approach 
of Sec. IIVI with standard puncture initial data. How- 
ever, we first need to inform the reader of a modification 
to the transition function in fl7| used for smoothing the 
far-zone and near-zone metric contributions in the out- 
ermost buffer zone. The parameters tq and w, which 
define where the transition begins and the width of the 
transition window respectively, have been adjusted. The 
analysis in [TtJ uses r = Xqw/5 and w — Xgw [see 
Eqs. (8.4) and (8.6) in 17|], which then put the mid- 
point (where the transition function has a value of 1/2) 
at r w 14Xgw- In the work here we choose values of 
r = 0.044Agw and w = 0.22AgWj which lead to a mid- 
point at r ps Xqw- This adjustment is justified because 
one would expect the transition region (from near zone to 
far-zone) to be more suitable closer to Xgw than IAXqw- 



1. Waveforms 

We first discuss the gravitational waveforms obtained 
from these data sets. In the figures that follow we de- 
note our results as Harmonic-match or ADMTT-PN; de- 
pending on whether we evolve initial data constructed by 
asymptotic matching of a second-order post-Newtonian 
4-metric expansion in Harmonic gauge to a perturbative 
black hole solution, or data from a post-Newtonian ex- 
pansion in ADMTT gauge. We do not expect to see the 
same deficiencies as seen in ADMTT-match. Figure [5] 
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FIG. 5: The (2,2) mode of ^4 extracted at a radius of r = 90M 
for standard punctures, Harmonic-match, and ADMTT-PN 
data sets. 

shows the real part of the (£ — 2,m — 2) waveform 



mode of ipi at an extraction radius of 90M for all three 
data sets. Harmonic-match produces gravitational waves 
starting from t = 0, more characteristic of an astrophys- 
ically realistic evolution. The diminished amplitude in 
the spurious radiation region for Harmonic-match is an 
additional positive outcome for this data set (although 
the ultimate goal is to completely eliminate this behav- 
ior from the waveform) . Figure [5] also shows the wave- 
form for ADMTT-PN and standard punctures. As with 
Harmonic-match, ADMTT-PN shows a non-zero wave- 
form signature at t = and comparable reduction of the 
spurious radiation amplitude. It is also noted that a com- 
parative' reduction of the spurious radiation amplitude 
has been seen in f44| as well. 



2. Eccentricity and Horizon Mass 

We can see from Fig. [5] that the waveform behav- 
ior of both Harmonic-match and ADMTT-PN is non- 
monotonic in amplitude. This differs from the standard 
chirp signal expected during the inspiral phase of the two 
black holes. The reason for this behavior is due to the 
black hole orbits being noticeably eccentric. The high 
eccentricity can be seen through plots of the coordinate 
separation over time (See Fig. [5]). For a sense of scale, 
the coordinate separation for standard puncture initial 
data is included to indicate the small level of deviation 
for an evolution to be accepted as having low eccentric- 
ity. Note, however, that unlike in the case of standard 
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FIG. 6: Coordinate separation for standard punctures, 
Harmonic-match, unmodified ADMTT-PN, and modified 
ADMTT-PN. 



punctures, the tangential momentum of the black holes 
has not been fine-tuned and is simply coming from PN 
expressions for circular orbits. For example, in the case 
of ADMTT-PN we use the momentum for circular PN 



8 



orbits given by 



PPN 



[Mil] 



6 



+e 



(441 -324r7-r7 2 )[Am] 
72 



(26) 



where M = mi + m,2, [i = mxmijM , rj = fi/M and f2 is 
defined as 



n 



i 

U 



64(^) 3 



(l + 2(i 7 )) 6 (m) 4 



(-0.625^ + ?y 2 

v Ml 



(27) 

Increased eccentricity, with respect to standard punc- 
tures, is again apparent in Fig. [7] for both the Harmonic- 
match and ADMTT-PN sets. The eccentricity calcula- 
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FIG. 7: Eccentricity estimation for standard punctures, 
Harmonic-match, unmodified ADMTT-PN, and modified 
ADMTT-PN. The range between 100M and 600M is the 
range to analyze all sets simultaneously while guaranteeing 
no errors due to incomplete information (Before 100M our 
eccentricity estimator which averages over an orbital time is 
unreliable, and after 600M most sets have begun merger so 
that one cannot define an eccentricity.). 

tion is performed as in (45[, and is calculated through a 
time-average over the orbital period. 

For the case of ADMTT-PN, we have added an ad- 
ditional parameter, C3, that can be used to adjust the 
initial tangential momenta of the black holes. The mod- 
ified tangential momentum, p is given by 



2 2^ 6 M 2 M M 

V =Ppn +e -z C 3 — 

8ri2 V r i2 



(28) 



Note that the term with C3 comes into the expression at 
order 0(v/c) 6 . Thus this additional term is of the or- 
der of one of the PN terms that were neglected in the 
expansion in Eq. (f26|). By adjusting C3 we can then de- 
crease the eccentricity. For C3 = —396.4844 we change 
our tangential momentum to 98.3% of the original ppn- 
This choice effectively minimizes the eccentricity. In 



Figs. [6] and [7] we can compare the results with and with- 
out this modification to Harmonic-match and standard 
punctures. 

A useful local definition for the mass of a black hole 
is the apparent horizon mass. An algorithm for finding 
apparent horizons without any symmetries is described 
in [461 ] . Over time, numerous variations on Gundlach's 
methodology have produced routines for many of the 
main computational groups involved in black hole sim- 
ulation. The version used here is known as AHmod and 
was developed by Norbert Lages {47} as an improvement 
to the CACTUS thorn AH finder developed by Miguel 
Alcubierre [48]. Figure [5] shows the apparent horizon 
mass of one of the (equal mass) black holes for all three 
data sets versus time. We see that ADMTT-PN and 
Harmonic-match show much more variation than stan- 
dard puncture data. We see that the two alternative 
methods lead to a "mass loss" over the evolution time 
where they eventually approach the apparent horizon 
mass obtained from evolving standard puncture data. 
This "mass loss" , albeit of low scale, is unphysical and 
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FIG. 8: Apparent horizon mass of one of the black holes in 
the inspiral phase under the standard punctures, Harmonic- 
match, and ADMTT-PN approach. 



we expect it is due to larger Hamiltonian constraint vio- 
lations in the approximate initial data sets. For example, 
in [44j the Hamiltonian constraint (but not the momen- 
tum constraint) is solved numerically for ADMTT-PN, 
which leads to a horizon mass that is closer to the result 
for standard puncture data. 

Figure [9] shows the apparent horizons masses of the 
final black hole after merger for evolutions with each 
data set. The masses obtained with the alternatives tend 
to deviate from punctures, with Harmonic-match yield- 
ing an increase and ADMTT-PN yielding a decrease in 
the common horizon mass respectively over time. This 
decrease is certainly unphysical and is likely related to 
larger constraint violations. 



9 



0.886 
0.8855 
^ 0.885 

s 

en 0.8845 
in 
a 

* 0.884 
O 

■= 0.8835 
O 
X 

0.883 
0.8825 

700 



Standard Punctures ■ 
Harmonic-match 
ADMTT-PN 




800 



900 



t(M) 



1000 



1100 



FIG. 9: Apparent horizon masses of the final black hole after 
merger under the standard punctures, Harmonic-match, and 
ADMTT-PN approach. 



3. Constraints 

The violation of the Hamiltonian constraint equation 
at the beginning of the evolution, as seen in Fig.fTOlis be- 
tween 2 and 3 orders of magnitude larger for Harmonic- 
match and ADMTT-PN than for standard punctures at 
t = 0. We expect that these violations are responsible 
for the observed drift in the apparent horizon masses. 
Over time, the evolution of ADMTT-PN relaxes the con- 
straints to the regime of standard punctures where both 
are almost indistinguishable. This result brings the ex- 
pectation that these initial data methods are stable over 
the life of the evolution, and a case for astrophysical rel- 
evance. 
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FIG. 10: The L -norm of the Hamiltonian constraint violation 
for the standard punctures, Harmonic-match, and ADMTT- 
PN approach. The y-axis is plotted in logio- The inset shows 
the violations at the beginning of the evolution. 

Figure QT] shows the L 2 -norm of the x-component of 
the Momentum constraint. It is seen that the initial 
discrepancies (compared to standard punctures) are of 



2 orders of magnitude for ADMTT-PN and roughly 5 
orders of magnitude for Harmonic-match. The y and z 
components are so similar that they would yield almost 
the same plots. Over time, all data sets relax to within 
an order of magnitude difference until they become al- 
most indistinguishable during the merger and ringdown 
phases. 
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FIG. 11: The L -norm of the x-component of the Momentum 
constraint violation, for the standard punctures, Harmonic- 
match, and ADMTT-PN approach. The y-axis is plotted in 
logio. The inset shows the violations at the beginning of the 
evolution. 



VII. DISCUSSION 

We have investigated three alternative approaches of 
constructing initial data for a binary black hole system. 
Two of these alternatives involve asymptotic matching 
of post-Newtonian expansions to perturbative black hole 
solutions of the metric, and one uses a resummed post- 
Newtonian expansion of the metric everywhere. We nu- 
merically evolve all three alternatives using the BSS- 
NOK system fll| together with the moving punctures 
approach 0, Qand compare them with evolutions that 
start from standard puncture initial data. The two alter- 
native data sets that were derived by asymptotic match- 
ing have not been evolved before. Since they contain 
singularities near the black hole centers, we have also 
developed a particular algorithm to fill a certain region 
inside each black hole with smooth data that approaches 
puncture data at the center. In this way we can avoid the 
use of black hole excision, so that we are able to use the 
same moving punctures approach to evolve all our data 
sets. 

We find that the simplest data set, using low-order 
asymptotic matching, does not lead to physically accu- 
rate evolutions. Insufficient expansion of the extrinsic 
curvature leads to black hole momenta that are far too 
small for quasicircular orbits, so that we observe an al- 
most head-on collision. Also, low-order expansion of the 
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metric yields gravitational waveforms lacking an initial 
signature. Thus, this initial data set is not an improve- 
ment over standard punctures. 

Higher-order expansion techniques resulted in the pro- 
duction of gravitational wave signals at the initial time 
as well as proper quasicircular inspiral orbits that subse- 
quently lead to merger and ringdown. Eccentricity mea- 
sures were higher for the Harmonic-match and ADMTT- 
PN sets than what is seen in standard punctures, when 
the puncture momentum parameters are fine-tuned to 
achieve low eccentricities. However, by tuning our 3PN 
order C3 parameter (see Eq. (l28|) ) we were able to bring 
the eccentricity of ADMTT-PN down as well. Both 
Harmonic-match and ADMTT-PN sets look promising in 
the sense that, unlike for standard puncture data, they 
have realistic gravitational wave signals built-in from the 
beginning of the evolution. In addition, they show a 
reduction of spurious radiation when we compare with 
standard puncture data. 

There are, however, still problems with these alterna- 
tive data sets that stem from the fact that they only 
approximately satisfy the constraint equations of gen- 
eral relativity. We have monitored the Hamiltonian and 
momentum constraint violations during the evolutions. 
We find that constraint violations of Harmonic-match 
and ADMTT-PN sets are above the usual level com- 
ing from numerical truncation errors observed when we 
evolve standard puncture data. In the case of ADMTT- 
PN these violations eventually decrease and reach a level 
similar to what is seen for standard puncture data evo- 
lutions. However, until this happens we see an unphys- 
ical drop in the apparent horizon masses. In the case 
of Harmonic-match, the constraint violations start at an 
even higher level and never quite fall to the level observed 
for standard puncture data evolutions, and again result 
in unphysical behavior of the apparent horizon masses. 

Future improvements need to be focused on these 
constraint violations. An obvious remedy would be to 
project the alternative data sets onto the solution mani- 
fold of general relativity by solving the constraint equa- 
tions, e.g. as in [l2j by using the York-Lichnerowicz con- 
formal decomposition (49|. If the constraints are satis- 
fied we expect the apparent horizon masses to be strictly 
non-decreasing over time. Evolutions, horizon masses, 
and eccentricity could then be compared to the standard 
punctures approach and conclusions could be made on 
more stringent grounds. 

Most astrophysical black holes are expected to be spin- 
ning. It is thus important to include spin in the initial 
data construction. As already mentioned in Sec. lIIIi stan- 
dard punctures contain spin parameters that enter the 
Bowen-York extrinsic curvature in Eq. ([5]), so that one 
can set up initial data with spins. Including spins in the 
alternative data sets discussed in this paper is possible 
as well. For any of the data sets that are produced from 
asymptotic matching one would have to redo the match- 
ing calculation as in [l7| and match a post-Newtonian 
4-metric with spin to a perturbed black hole metric with 



spin. This calculation, while certainly long and tedious, 
is in principle possible. The situation is much simpler 
for the post-Newtonian data set in ADMTT gauge dis- 
cussed in Sec. IIVI To leading order, the spin only enters 



the Bowen-York piece A l ^ Y in Eq. (|15l) in the same way 
as in the case of standard punctures. It would thus be 
easily possible to add these missing spin terms. This 
alone would already result in spinning black holes. How- 
ever, in order to obtain the correct waveform at the ini- 
tial time we have to ensure that i n Eqs. Q13[) and 
(fTS"]) is computed for post-Newtonian particle trajecto- 
ries where spin is included in the equation of motion. 
Since post-Newtonian trajectories for spinning particles 
are well known, and since is known for arbitrary 
trajectories [HI], the inclusion of spin in this data set is 
not a complicated problem. 
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Appendix A: Drawbacks of Asymptotic Matching 



In [151 - 1171 the ansatz, 
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X'tf) = £ (^y 72 (X a )j (x*) + 0(v/cf (Al) 

is made for the coordinate transformation that trans- 
forms the coordinates X a of the perturbed black hole 
solution in inner zone 1 to the coordinates x" that are 
used in the near zone for the post-Newtonian metric. 

Note that (2f-) j/2 



17 it is 



0(v/ C y. In Sec. V of 
shown that the zeroth order piece {X a ) of this coordi- 
nate transformation corresponds to a simple translation 
and the first order piece is 



{X a ) 1 = (i?> C( ) 1 ^ + (C Q ) 1 , 



(A2) 



where (Fp a ) 1 is anti-symmetric. The (Fp a ) 1 and (C a ) 1 
correspond to constants of integration that appear when 
matching at 0(v/c). As shown in Sec. V.E. of [T3|, the 
(-F/3a)i can only be determined by matching the two 4- 
metrics up to 0(v/c) 3 . To be precise, (Fp a ) 1 comes 
from matching the and goo components at 0(v/c) 3 . 
This pattern continues in the sense that at each order 
of matching one finds free constants of integration that 



11 



can only be fixed once matching has been performed two 
orders higher in v/c. 

In the case of the lower order matching procedure used 
in (l5l . [l6j , matching of the and goo components of the 
two 4-metrics has only been performed up to 0(v/c) 2 . 
This means that the piece {Fp a ) 1 could not be computed. 
In fact, it was merely set to a value that simplified fur- 
ther calculations. This has important consequences. As 
shown in Sec. V.E. of [17| the only non-zero components 
of {Fp a ) x are (F 02 ) 1 = - (^20)1 = -y/m 2 /m. This yields 

a term (^) 1/2 (F^ t = {ff'tin Y = X\x?), 
which corresponds to a boost. This boost is needed to 
transform the metric of the perturbed black hole solu- 



tion in inner zone 1 (from a frame where it is at rest to a 
frame where it is moving in t he ^ -direction) . In the case 
of lower order matching as in p~HL [TgI] this term is absent, 
so that the metric used near the black hole is missing this 
boost. The resulting extrinsic curvature (which is always 
smaller than the 3- metric by a factor of v/c) is thus cor- 
rect only up to 0(v/c) 2 and has errors of 0(v/c) 3 , where 
the missing boost would enter. This leads to a black hole 
without sufficient momentum. 

In general, this implies that if we want to obtain a 
4-metric that is correct up to 0(v/c) n with the match- 
ing procedure in [lB - fi~7j we need to find the matching 
coordinate transformation up to order 0(v / c) n+2 . 
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